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ABSTRACT 


The non-detection of a coma surrounding 1I/‘Oumuamua, the first discovered interstellar object 
(ISO), has prompted a variety of hypotheses to explain its nongravitational acceleration. Given that 
forthcoming surveys are poised to identify analogues of this enigmatic object, it is prudent to devise 
alternative approaches to characterization. In this study, we posit X-ray spectroscopy as a surprisingly 
effective probe of volatile ISO compositions. Heavily ionized metals in the solar wind interact with 
outgassed neutrals and emit high-energy photons in a process known as charge exchange, and charge 
exchange induced X-rays from comets and planetary bodies have been observed extensively in our 
Solar System. We develop a model to predict the X-ray flux of an ISO based on its chemical inventory 
and ephemeris. We find that while standard cometary constituents, such as H20, CO», CO, and dust 
are best probed via optical or infrared observations, we predict strong X-ray emission generated by 
charge exchange with extended comae of Hz and Nz — species which lack strong infrared fluorescence 
transitions. We find that XMM-Newton would have been sensitive to charge exchange emission from 
1I/*Oumuamua during the object’s close approach to Earth, and that constraints on composition may 
have been feasible. We argue for follow-up X-ray observations of newly discovered ISOs with close-in 
perihelia. Compositional constraints on the general ISO population could reconcile the apparently self- 
conflicting nature of 11/‘Oumuamua, and provide insight into the earliest stages of planet formation 


in extrasolar systems. 
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1. INTRODUCTION 


Knowledge of the compositions of minor bodies has 
considerably influenced our understanding of the early 
Solar System. While investigations of their dynamical 
histories are limited by the Solar System's chaotic na- 
ture (Wisdom 1980; Laskar 1989; Batygin & Laughlin 
2008; Laskar & Gastineau 2009), compositional mea- 
surements of minor bodies — especially of their interiors 
— have the potential to reveal their formation environ- 
ments (Oberg et al. 2011). 

Remarkably, the Solar System likely ejected tens of 
earth masses of volatile rich material during the early 
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stages of its formation and evolution (Hahn & Malhotra 
1999; Gomes et al. 2004; Tsiganis et al. 2005; Morbidelli 
et al. 2005; Levison et al. 2008; Raymond et al. 2018, 
2020). Therefore, it is reasonable to expect that extraso- 
lar planetary systems contribute to a galactic population 
of interstellar comets, and that some will encounter the 
Solar System on hyperbolic trajectories (Moro-Martin 
et al. 2009; Engelhardt et al. 2017; Cook et al. 2016). 
Naturally, it came as a surprise that 1I/‘Oumuamua, 
the first discovered ISO, exhibited none of the typical 
properties of Solar System comets (for recent reviews, 
see Moro-Martin 2022; Jewitt & Seligman 2022). It had 
no visible coma (Meech et al. 2017; Jewitt et al. 2017; 
Trilling et al. 2018; Bolin et al. 2018), an extreme shape 
(e.g., Knight et al. 2017; Mashchenko 2019), a reddened 
reflection spectrum (Masiero 2017), a young dynamical 
age (e.g., Mamajek 2017; Hsieh et al. 2021) and non-zero 
nongravitational acceleration (Micheli et al. 2018). 


The second discovered ISO, 2I/Borisov, was unam- 
biguously a comet (Jewitt & Luu 2019), perhaps more 
inline with expectations. Its cometary activity was read- 
ily measured and typical carbon and nitrogen bearing 
species were detected (Opitom et al. 2019; Kareta et al. 
2020; Lin et al. 2020; Bannister et al. 2020; Xing et al. 
2020; Aravind et al. 2021). ALMA and Hubble observa- 
tions revealed that it was enriched in CO relative to H5O 
(Bodewits et al. 2020; Cordiner et al. 2020). This find- 
ing suggests that 2I/Borisov formed exterior to the CO 
snowline in the protoplanetary disk (Price et al. 2021; 
Seligman et al. 2022) of a very young (Lisse et al. 2022) 
or a very carbon-enriched system (cf. Bodewits et al. 
2020). 

‘Oumuamua’s former provenance remains in question, 
and dynamically backtracing either ISO to their respec- 
tive host systems is difficult (Mamajek 2017; Bailer- 
Jones et al. 2018). 1I/‘Oumuamua’s peculiar prop- 
erties and elusiveness to follow-up spectroscopy have 
prompted a variety of theories regarding its origin. If 
the measured nongravitational acceleration was driven 
by cometary outgassing, this would be energetically con- 
sistent with a composition of Hə (Füglistaler & Pfen- 
niger 2018; Seligman & Laughlin 2020; Levine & Laugh- 
lin 2021), Ng (Jackson & Desch 2021; Desch & Jackson 
2021), or CO (Seligman et al. 2021). Confirming an Ha 
or Ng composition would have immediately revealed a 
new mechanism of generating minor bodies: formation 
in the cores of molecular clouds (Seligman & Laughlin 
2020), or ejection of ice from the surfaces of Pluto ana- 
logues due to impact events (Desch & Jackson 2021). 
Alternatively, if the acceleration was caused by radia- 
tion pressure (Micheli et al. 2018), then 1I/‘Oumuamua 
must have been extremely porous (Moro-Martín 2019; 
Sekanina 2019; Flekkøy et al. 2019; Luu et al. 2020) or 
thin (Bialy & Loeb 2018). While H20, CO2, and CO 
species in cometary comae are easily measured via in- 
frared spectroscopy, the possible presence of difficult to 
detect homonuclear Hz or N2 prompts a search for al- 
ternative methods for characterizing interstellar comets. 

This study presents such a new approach: X-ray ob- 
servations of solar wind (SW) charge exchange (CX). 
X-ray observations of an interstellar object should im- 
mediately reveal the presence of an outgassed coma. 
CX involves the transfer of electrons from a cool neu- 
tral medium to heavily ionized metals in the SW. Fol- 
lowing the unexpected detection of X-rays emanating 
from C/1996 B2 (Hyakutake) (Lisse et al. 1996; Cravens 
1997), CX has been observed with a multitude of other 
comets (Dennerl et al. 1997; Cravens 2000, 2002; Lisse 
et al. 1999; Lisse et al. 2004; Bodewits et al. 2007). CX 
conveys information about the composition and speed 


of the solar wind (Beiersdorfer et al. 2001; Bodewits 
et al. 2007; Gu et al. 2016). Studies have also linked 
CX emission to physical and chemical properties of the 
minor body and its coma through the morphology of 
the X-ray emission profile (Wegmann et al. 2004) and 
emission line ratios (Mullen et al. 2017). 

Indeed, there are encouraging prospects for discov- 
ering additional ISOs in the coming decade. Popu- 
lation synthesis models of the galactic population of 
ISOs have demonstrated that the Rubin Observatory 
Legacy Survey of Space and Time (LSST) will detect 
1—2 11/‘Oumuamua analogues every year (Hoover et al. 
2022). 

This paper is organized as follows. In 82, we dis- 
cuss the advantages and limitations of X-ray observa- 
tions for characterizing minor bodies, and outline a sim- 
ple model of charge exchange emission. We apply our 
model to 11/‘Oumuamua and 2I/Borisov in 83, and com- 
pare our predictions to the known correlation between 
X-ray and visible-band flux. Our model suggests that 
if 11/‘Oumuamua did indeed outgas a coma, that CX 
would have been detectable with current-generation X- 
ray facilities during the period that it approached Earth. 
§4 establishes expectations for detecting CX with new 
ISOs discovered by the Rubin Observatory. Finally, our 
conclusions are summarized in 85. We refer the reader 
to Appendix A for a brief review of CX fundamentals, 
and to Appendix B for the details of our CX model. 


2. MOTIVATION & MODEL 


This section establishes use cases for X-ray observa- 
tions of interstellar objects, and sets forth an analytic 
model of X-ray emission from outgassed comae. It also 
specifies the traits of ISOs (i.e., bulk composition and 
orbital trajectory) from which we expect robust detec- 
tion of X-ray emission with present-day facilities. 


2.1. Information Gained from X-ray Observations 


As a primer, it is useful to review the prototypical case 
of CX from an outgassed coma: C/1996 B2 (Hyaku- 
take), a long-period comet which came within 0.23 au 
of the Sun at perihelion. Infrared observations yielded 
unambiguous detections of H20, CO, CH4, and C2H6 
(Mumma et al. 1996; Dello Russo et al. 2002), Con- 
temporaneously, significant X-ray emission was detected 
out to ~ 100,000 km (6.7x107~* au) from its nucleus 
by Lisse et al. (1996). CX was suggested as the re- 
sponsible process by Cravens (1997), and reaffirmed by 
Wegmann et al. (1998). Other mechanisms (e.g., ther- 
mal bremsstrahlung, scattering, and fluorescence) may 
be ruled out (Krasnopolsky 1997; Lisse et al. 1999; Lisse 
et al. 2004). Since this detection, sophisticated models 


of CX with cometary comae have been developed, such 
as the hydrodynamical simulations by Wegmann et al. 
(2004), and CX has been observed with comae of nu- 
merous other comets (e.g. Ewing et al. 2013; Bodewits 
et al. 2007). 

X-rays have been conventionally of limited value for 
compositional studies of minor bodies. For one, CX de- 
tections arise only from significantly extended comae. 
Moreover, the most immediate insights conveyed are 
the morphology of the coma and the properties of the 
solar wind (Lisse et al. 2004). Infrared/radio spec- 
troscopy remains the de facto technique for detecting 
molecular species which exhibit solar-pumped fluores- 
cence transitions (Crovisier & Encrenaz 1983; Bockelée- 
Morvan et al. 2004). Such species (reviewed by e.g. 
Biver & Bockelée-Morvan 2016; Bockelée-Morvan & 
Biver 2017) include H30, CO, CO», CHa, C3H5, CoHe, 
and HCN. Complementary to infrared observations, UV 
spectroscopy can measure the production of some addi- 
tional species (Biver et al. 2022) such as noble gases and 
homonuclear diatomic molecules. However, except for a 
small number of tenuous reports of their sublimation 
(e.g., Ar in Hale-Bopp, Stern et al. 2000) noble gases 
remained elusive in comets prior to in situ detections in 
67P (Rubin et al. 2018) by the European Space Agency’s 
Rosetta mission. Rosetta also confirmed the presence O» 
and N» in comets (Rubin et al. 2015; Bieler et al. 2015). 

This study explores a novel use case for X-ray obser- 
vations: to reveal an outgassed coma that is invisible in 
the infrared, and which may arise from the perihelion 
passage of an interstellar object. While the composition 
of 1I/‘Oumuamua remains debated, a number of candi- 
date species (Table 1) are monotomic or homonuclear di- 
atomic, and will not fluoresce in the infrared. However, 
CX will occur at a neutral medium, even one that does 
not produce observable thermal emission, fluorescence, 
or scattering. Therefore, the outgassed coma of an inter- 
stellar object may be X-ray bright even in the absence of 
infrared and optical detections. In particular, X-ray ob- 
servations offer a promising avenue for remote sensing an 
Hə or Nə coma. Both of these species are thermodynam- 
ically consistent with 1I/*Oumuamua's nongravitational 
acceleration under outgassing, and comprise the focus 
of this study. Per Seligman et al. (2021), we benchmark 
results to a CO composition (c.f. Trilling et al. 2018); 
although this molecule is more readily detected in the 
infrared. 


2.2. Scaling Relationships 


In Appendix B, we develop a model that links X-ray 
flux to a minor body's physical properties. Specifically, 
it depends on the following key parameters: the neutral 
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Table 1. Cryogenic chemical properties for volatiles which 
may be constituents of ISOs. For each species, we report 
the temperature of sublimation (Tsun), ideal adiabatic index 
(y), enthalpy of sublimation (AH), and sound speed (cs) per 
Equation B13. We also report the total energy input to a 
single particle (H) per Equation B12. Data from Shakeel 
et al. (2018), except for CO which is obtained from Stephen- 
son & Malanowski (1987) and the NIST database. 


Species Tsun ¥ AH Cs H 
(K) (kJ mole?) (kms !) (kJ x 107°) 
Hə 6 7/5 1 0.19 0.18 
Ne 9 5/3 1.9 0.11 0.34 
N2 25 7/5 7.34 0.14 1.27 
CO 58 7/5 7.6 0.15 1.37 
Ar 30 5/3 7.79 0.15 1.36 
O2 30 7/5 9.26 0.10 1.60 
Kr 40 5/3 11.53 0.12 2.01 
Xe 55 5/3 15.79 0.12 2.75 
CO. 82 8/6 28.84 0.14 4.94 
H:O 155 8/6 54.46 0.31 9.33 


production rate (Qgas), mean radius of the nucleus (a), 
perihelion distance (rperi), sublimation + kinetic energy 
of the relevant volatiles (H), geocentric distance (Ae), 
solar wind number density (ngw), CX interaction area 
of the coma (A), and X-ray lumonisty (Lx) and flux 
(Fx). The pertinent relationships are as follows: 


Q gas xa je Ma (1) 
A x Q2as (to leading order) (2) 
nsw xX pne (3) 
Lx x nswA (4) 
Fx x A;?Lx « atr SH? Ae (5) 


In our model, constants of proportionality are cal- 
ibrated to measurements of C/1996 B2 (Hyakutake) 
(Lisse et al. (1996), Appendix B). The above scaling 
relationships demonstrate that the X-ray brightness of 
a minor body can be significantly enhanced by a close-in 
perihelion passage and a particularly volatile composi- 
tion. Figure 1 shows a schematic of such an ISO en- 
counter, where the scaling relationships and underlying 
physical processes are depicted in the individual panels. 


2.3. Detection Thresholds 


In this section we guage the feasibility of detecting 
CX with an outgassed coma, including the requisite X- 
ray flux and exposure durations. Our predictions are for 
the EPIC pn detector onboard XMM-Newton, which has 
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Figure 1. Schematic of solar wind charge exchange with the outgassed neutrals in the coma of an ISO. Panel 1: An instantaneous 
view of the inner solar system and the trajectory of an ISO, currently at perihelion (rperi). Panel 2: The flux of solar radiation 
(Fo) heats the volatiles on the ISO's surface, which causes sublimation and outgassing of neutral volatiles at a rate Qeas. 
Panel 3: Solar wind ions (with number density ngw) interact with the outgassed coma through charge exchange. Electrons are 
transferred from the neutral coma to the bare or H-like solar wind ions, and the subsequent cascade emits X-rays. The coma’s 
interaction cross-section (A) scales with Q2,.. As a result, the X-ray luminosity follows Lx œ r5°, (Equation 5, 82.2). 


the advantage of a larger effective area (approximately 
A = 900 cm? at 0.5 kev) compared to the other space- 
based X-ray detectors. To estimate the detectability of 
a comet or ISO, one needs to account for the noise con- 
tribution from the background, which is chiefly made of 
background AGNs and the diffuse hot plasma emission 
from the general interstellar and circumgalactic media 
of the Galaxy, the Local Bubble, and solar wind charge 
exchange emission within the heliosphere (e.g., McCam- 
mon et al. 2002; Koutroumpa 2012). Based on archival 
blank-sky EPIC pn observations, we estimate that the 
total background intensity is © = 1.4 x 107? counts s^! 
arcmin~? in the 0.5 — 0.7 keV band, for example, which 
encloses the typically strong Ovir and Ovin Ka lines. 
The total background intensity also includes the contri- 
bution from non-X-ray (cosmic-ray induced instrument) 
events. Let ®x denote the X-ray flux at the detector due 
to CX with the target (units of counts s^! cm~?), and 
let + denote the exposure duration. Furthermore, we 
assume that nearly all of the signal is contained within 
a solid angle 2 for simplicity. Then the signal-to-noise 
(S/N) of the CX signature is 


S/N = 6x AV/r/ / bx A+ EQ. (6) 


Alternatively, this equation can also be rearranged into 
an expression for 7 for a desired S/N threshold. 


peri 


One caveat is that a continuous exposure may be inef- 
fective for fast-moving targets. Photon events may, how- 
ever, be spatially calibrated across discrete slews that 
keep the object on the focal plane array (FPA). This 
enables multiple on-target pointings and reconstruction 
of events in the cometocentric frame based on the atti- 
tude of the telescope. The maximum exposure duration 
for any one pointing is then limited by the field of view 
of the telescope and the angular velocity of the target. 
For reference, 11/^Oumuamua's angular velocity reached 
1.4 x 1074 deg. s-! at closest approach to Earth. The 
crossing time for XMM-Newton’s 30' x 30' field of view 
would have been up to ~ 3.5 ks. 

Space-based X-ray observatories have minimum so- 
lar elongation angles (7%) for reasons pertaining to 
detector safety (e.g., avoiding burn-out incidents simi- 
lar to ROSAT's), thermal stability, and power supply. 
For Chandra, min = 46.4°, and for XMM-Newton, 
ws" = 70°. These constraints permit observations of 
small bodies at R > 0.72 au and R > 0.94 au, respec- 
tively. These lower bounds can only be attained un- 
der optimal alignment between Earth and the object. 
Future missions will likely have similar safety exclusion 
zones to protect equipment at the FPA. 

Figure 2 depicts the exposure times necessary to 
achieve S/N = 5 as a function of Qgas and rperi, Specifi- 
cally for the XMM-Newton pn detector. For simplicity, 
the geocentric distance is assumed A, = 1 au. The 
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Figure 2. Exposure time (ks) required for a 5e detection 
of CX X-ray emission (Equation B15), as a function of out- 
gassing rate Qgas and perihelion distance rperi. Geocentric 
distance is fixed at Ae = 1 au. The background rate and ef- 
fective area correspond to that of EPIC pn on XMM-Newton. 
For reference, markers are plotted at locations representa- 
tive of real minor bodies: light-blue (‘H’) corresponds to 
Qeas = 1.5 x 107? s^! observed at 1.0 au — similar to C/1996 
B2 (Hyakutake); tan (‘B’) corresponds to Qgas = 10?" s^! 
observed at 2.0 au — similar to 2I/Borisov; coral (‘O’) cor- 
responds to Qaa, = 2.8 x 10°85 s^! observed at 0.25 au 
and Qgas = 1.8 x 10°” s^! observed at 1.0 au — estimates 
corresponding to 1I/*Oumuamua under a CO composition. 
The solar elongation angles of Chandra (i$" = 46.4?) and 
XMM-Newton (3 = 70°) impose lower bounds (sin 8") 
on observable rj, and their respective restricted zones are 
hatched in the figure. The lower bound is attained only in 
the ideal orbital case and observable heliocentric distances 
are typically significantly larger due to spacecraft limitations. 


X-ray luminosity of CX is given by Equation B15, and 
the angular extent of the emission is set by b (Equa- 
tion B10). Markers in Figure 2 are located at outgassing 
rates and perihelia of real minor bodies (the outgassing 
rate of 11/‘Oumuamua was not directly measured, but 
is estimated in §3). 

This figure immediately conveys the feasibility of de- 
tecting CX emission from outgassed comae with cur- 
rently operational facilities. For modest exposure times 
of order 10 ks, CX should be detectable when Q > 10?7 
s^! over a range of rperi- Qualitatively, this finding is 
consistent with previous high-significance ROSAT detec- 
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tions of CX with six comets (Dennerl et al. 1997), whose 
ougassing rates were between Qgas = 2 x 10?” — 9 x 109 
s71. The heliocentric distances in this sample were be- 
tween R = 0.99 — 1.96 au, and geocentric distances be- 
tween A, = 0.12 — 1.61 au (these comets are discussed 
further in §3.3). While the minor bodies need not be ob- 
served at perihelion, our model indicates that Fx scales 
more strongly with heliocentric distance (R) than with 
Ae (Lisse et al. 1999) per Equations B11 & B14. As 
discussed in the next section, 1I/‘Oumuamua’s proxim- 
ity to the Earth and Sun would have made it amenable 
to X-ray observations. On the other hand, 2I/Borisov 
would have been a more difficult case because it only 
came within a minimum 2 au of the Sun. 


3. DETECTABILITY OF X-RAYS FROM 
2I/BORISOV AND 1I/'OUMUAMUA 


We proceed to estimate the X-ray luminosity for the 
only two identified ISOs that traversed our Solar System: 
1I/‘Oumuamua and 2I/Borisov. We first estimate Lx 
using the analytic framework presented in the previous 
section. We then make an independent estimate based 
on known trends between optical and X-ray luminosities 
for comets. 

We briefly review pertinent measurements of both 
ISOs. From its photometric lightcurve structure, 
1I/*Oumuamua is estimated to have had an oblate 6:6:1 
ellipsoidal geometry (Mashchenko 2019). Its albedo (p) 
was not constrained, so its true dimensions are un- 
certain. Meech et al. (2017) and Bolin et al. (2018) 
measured an effective radius of ~ 0.1 km assuming 
an albedo of 0.04. The light curve fit by Mashchenko 
(2019) yielded dimensions 115 x 111 x 19 m with p = 
0.1. For our order-of-magnitude calculations, we treat 
1I/‘Oumuamua as spherical with radius of a = 100 m. 
However, we emphasize that the true dimensions are 
albedo-dependent, and the lengths scale as 1/,/p. At 
perihelion of ry; = 0.25 au, 11/‘Oumuamua’s distance 
to Earth was A, = 1.23 au with an apparent magni- 
tude V = 21. The object reached an apparent magni- 
tude of V = 19.7 at its closest approach to Earth with 
A. = 0.16 au. Infrared Spitzer observations produced 
non-detections when 1I/'Oumuamua was outbound at 
~ 2 au, providing upper limits on the production rates 
of CO or CO» (Trilling et al. 2018). However, sporadic 
outgassing of CO might reconcile the nongravitational 
acceleration and Spitzer non-detection (Seligman et al. 
2021). 

The nuclear radius of comet-like 21/Borisov was con- 
strained to 0.2km « a « 0.5km based on its bright- 
ness profile and nongravitational acceleration assuming 
a range of plausible bulk densities (Jewitt & Luu 2019; 
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Bolin et al. 2020). We adopt a nominal radius of a = 300 
m. Unlike 1I, gas emission activity typical of comets was 
definitely detected from this body. When contempora- 
neous measurements of H2O and CO production rates 
were obtained after perihelion (rperi = 2 au), CO dom- 
inated the outflow. Specifically Bodewits et al. (2020) 
measured Q(CO) = 7.5 x 1076 s^! at 3.3 days after 
perihelion and Q(CO) = 1.07 x 10?” s7! at 22.4 days 
after perihelion. Also, Cordiner et al. (2020) measured 
Q(CO) = 4.4 x 1076 s^! in observations conducted 7 — 8 
days after perihelion. While water production exceeded 
Q(H30) = 10?" s~! prior to perihelion, it dropped below 
6 x 1076 s71 for observations conducted after perihelion 
(Bodewits et al. 2020). Jewitt et al. (2020) measured 
V = 16.6 at 26 days after perihelion, when 2I/Borisov's 
distances were R — 2.1 au and A, — 1.9 au. 

In the following, Qgas denotes the total outgassing rate 
across all neutral volatiles. In certain instances, a dom- 
inant species X was determined, in which case its out- 
gassing rate is denoted Q(X). Our predictions for Lx 
neglect subtle dependencies on properties of the solar 
wind and coma. In the optically thick regime, all solar 
wind minor ions will undergo charge exchange. However, 
the X-ray lines emitted depend on the neutral compo- 
sition (Mullen et al. 2016). Also, CX cross sections de- 
pend strongly on the solar wind velocity (Kharchenko & 
Dalgarno 2001; Beiersdorfer et al. 2001; Bodewits et al. 
2007). 


3.1. Analytic Predictions for 2I/Borisov 


For 2I/Borisov, we adopt the measured production 
rate Q(CO) = 1.07 x 10?” s-! at R = 2.07 au (Bode- 
wits et al. 2020). We note that Equation B11 would 
predict Qeas = 5.9 x 107° s^! near perihelion, and mod- 
ifying the scaling relation to take into account a CO 
composition (i.e., an H appropriate for CO ice) yields 
Qeas = 4.0 x 1027 st, which are both within an order- 
of-magnitude of the observed value. Accounting for 
the reduced solar wind density at 2 au, we calculate 
Lx = 7.5 x 10! erg s^! (see Table 2). A 900 cm? ef- 
fective area (characteristic of EPIC pn) at Earth would 
detect 6.3 x 107? counts s^ !, which would require a con- 
siderable 400 ks exposure for a 5c detection. Indeed, 
most previous CX detections in comets have involved 
much brighter X-ray luminosities Lx > 10!^ erg s7! 
(Lisse et al. 1999; Lisse et al. 2004). 

Since there are observations of 2I/Borisov at multi- 
ple points along its trajectory, we may perform a more 
detailed comparison between our model and measure- 
ments of outgassing over time. We plot the observed 
time-evolution of Q(H5O) and Q(CO) in Figure 4 with 
our production rate model (Equation B11). Our model, 


which was calibrated to observations of C/1996 B2 
(Hyakutake), predicts water production rates within a 
factor a few of those observed. The discrepancies may 
be due to composition inhomogeneities in the nuclues 
and/or the sporadic nature of outgassing. On the other 
hand, our high model values show that Borisov could 
not be made of pure CO, but is more in line with the 
~ 10— 2096 vs. water maximal abundances of CO found 
in solar system comets (Bockelée-Morvan & Biver 2017). 


3.2. Analytic Predictions for 11/*Oumuamua 


As a baseline estimate for 1I1/‘Oumuamua, we assume 
a bulk composition of CO, which is energetically com- 
patible with the nongravitational acceleration and re- 
vised Spitzer limits on the production rate given the 
low activity state outbound at 2 au (Trilling et al. 2018; 
Seligman et al. 2021). Equation B11 with H correspond- 
ing to a CO composition suggests Qgas = 2.8 x 10?? 
s^! at perihelion (the leftmost marker in Figure 2), and 
Qeas = 4.4 x 1076 s^! at 2 au, in agreement with the 
Seligman et al. (2021) estimate, as expected. 

With X-ray flux received at Earth depending strongly 
on the target's distance from the Sun, 1I/‘Oumuamua 
would have only become accessible to XMM-Newton 
at its minimum observable heliocentric distance R = 
0.943 au. This geometry is depicted in Figure 3. It 
would have been prudent to observe 1I/*Oumuamua in 
X-ray at this time, which roughly corresponded with 
1I/*Oumuamua's closest approach to Earth. 

At R — 0.943 au and A, — 0.276 au, we estimate 
Que = 2.0 x 107 s71, Lx = 1210" ees, anda 
flux at Earth of Fx = 5.4 x 10714 erg s^! cm™?. As- 
suming an Nə (Desch & Jackson 2021) or Hə (Seligman 
& Laughlin 2020) composition yields even higher Fx, 
as reported in Table 2. For any of these three highly 
volatile compositions, a signal would have been unam- 
biguously detected in a « 1 ks X-ray follow-up exposure, 
producing a unique, easily obtained, positive outgassing 
detection. 


3.3. Observations of CX with Comets 


Detections of CX from Solar System comets help 
gauge conditions under which we may detect CX from 
ISOs. The detection of X-rays from C/1996 B2 (Hyaku- 
take) (Lisse et al. 1996) prompted Dennerl et al. (1997) 
to search archival ROSAT data for other instances of CX 
with comet comae. Confirmed cases include: C/1990 
Kl (Levy), C/1990 Ni (Tsuchiya-Kiuchi), 45P (Honda- 
Mrkos-Pajduskova; HMP), and C/1991 A2 (Arai). The 
faintest target was 45P, with Fx = 1.5 x 10713 erg s^! 
cm”, which was detected at high confidence by ROSAT 
with a 191 s exposure. 


A 


‘Oumuamua at perihelion 


Figure 3. Trajectory and X-ray observability of 
1I/‘Oumuamua, prior to its discovery. Plotted lines connect 
the instantaneous positions of Earth and 11/‘Oumuamua. 
Blue lines denote a solar elongation angle meeting X-ray 
observing requirements, yo > v2", where yon = 70° for 
XMM-Newton. The smallest heliocentric distance allowed by 
this constraint is denoted R* = 0.943 au. At this time, the 
distance between Earth and 1I/*Oumuamua was A? = 0.276 
au. 


'This comet sample serves as an independent diagnos- 
tic for selecting ISOs with detectable CX emission. It 
is established that a correlation exists between comet 
X-ray (Lx) and optical (Lop) luminosities (e.g., Lisse 
et al. 1999; Lisse et al. 2004). Lx plateaus for the brigh- 
est comets, but follows linear trend for modest luminosi- 
ties Lx < 1015 erg s^! where Lx ~ 1074Lopt; however, 
there is dependence on the dust-to-gas ratio [D/G] that 
induces 1 — 2 dex of scatter. Using the trend presented 
by Dennerl et al. (1997), we estimate 2I/Borisov's X-ray 
luminosity as Lx = 6.1 x 101? erg s^! near perihelion, 
based on its V = 16.6. For 11/‘Oumuamua, at close ap- 
proach to Earth and V — 19.7, we find Lx — 2.6 x 10? 
erg s !. (Note that this estimate likely underestimates 
X-ray emission from 1I/*'Oumuamua, and serves only as 
a foil to other estimates presented in this section.) 

It is worth mentioning that dusty comae (e.g., 
C/Hale-Bopp 1995 O1 or 17P/Holmes, Lisse et al. 1997, 
2013) may destroy solar wind ions without concomi- 
tant CX X-ray production (Dennerl et al. 1997; Lisse 
et al. 2004). This effect may prove an obstacle to X- 
ray observations of dusty ISOs. However, the typical 
dust content of ISOs is unclear. Jewitt et al. (2017) 
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and Meech et al. (2017) placed stringent upper bounds 
on 1I/*Oumuamua's dust production rate, at 2 x 107+ 
kg s^! and 1.7 x 107? kg s^! respectively. Estimates 
of 2I/Borisov's reach as high as 35 kg s^! (Kim et al. 
2020). 


4. FUTURE ISO CX X-RAY EXPECTATIONS 


For the prospects of detecting CX X-ray emission with 
a future ISO, it is important to establish: (1) a set of 
criteria that provides high-odds of detecting CX X-ray 
emission; and (2) an estimated fraction of ISOs discov- 
ered that will exhibit a detectable X-ray flux. To ac- 
complish these goals, we conduct a Monte Carlo analy- 
sis assuming Rubin Observatory sky parameters (as the 
largest and deepest continual all-sky survey of the next 
decade, e.g., Li et al. 2022) which incorporates baseline 
predictions for the distribution of perihelia and geocen- 
tric distances of ISO trajectories, as well as their size 
distribution and possible compositions. 


4.1. Trajectories 


Our Monte Carlo analysis involves randomly drawing 
orbital configurations for ISOs from previous population 
synthesis studies. We predict X-ray flux at the minimum 
heliocentric distance that satisfies the solar elongation 
angle constraint of XMM-Newton. This alignment oc- 
curs at time £*, defined as: 


t* =argmin ||R(t)|| 
t 


1( (Q(0- R()-R() o 
dra ( IR RORO] 2 p 


where R and R’ are vectors representing the posi- 
tions of the ISO and Earth, respectively, with the Sun 
at the origin. The ISO’s hyperbolic trajectory may 
be expressed in a Cartesian coordinate system where 
R —(R, Ro, 0}: 


Rg = e2—1 
x y(Rı + eg/(e — 1)? — (a/(e - 1)? (8) 


with eccentricity e > 1 and perihelion rperi sampled 
from distributions provided by Hoover et al. (2022). 
Their study coupled dynamical simulations with the lo- 
cal stellar velocity dispersion to constrain the popula- 
tion of 11/‘Oumuamua-like ISOs that the Rubin Obser- 
vatory will find. They computed the probability density 
of perihelia, which indicates that 10%, 25%, and 50% of 
detected ISOs will attain rperi < 0.35, 0.54, and 0.83 
au, respectively. Under this distribution, it is clear that 
2I/Borisov had an unlikely large perihelion at 2 au (and 


mo € co OH 
ge es ee ee L L 
80 60 40 20 0 20 40 60 


Time Since Perihelion (days) 


Figure 4. Time-evolution of production rates for 2I/Borisov near perihelion (for CO, H2O, and OH). Figure adapted from 
Seligman et al. (2022). References for data are provided in their Table 2. The measurements demonstrate that CO production was 
substantially greater than H2O production after perihelion — the only time when both species were observed contemporaneously. 
Our production rate model (Equation B11) is plotted as dashed lines, holding 2I/Borisov’s radius as a constant a = 0.3 km, 
and adopting an appropriate H for H2O and CO. We obtained 2I/Borisov's heliocentric distance with a rebound (Rein & Liu 
2012) simulation using the Mercurius integrator (Rein et al. 2019). For reference, the trajectory of 21/Borisov through the Solar 
System is shown at right, with perihelion marked. The solid blue portion corresponds to the timespan plotted in the left panel. 


Table 2. Estimated X-ray luminosity (Lx) and flux (Fx) of 2I/Borisov and 1I/‘Oumuamua using our model. For 
11/‘Oumuamua, quoted outgassing rates (Qgas) depend on the object's assumed bulk composition, and heliocentric (R) and 
geocentric (Ae) distances. 11/‘Oumuamua’s pericenter was R = rperi = 0.25 au; however, the solar elongation constraint for 
XMM-Newton was first satisfied at R = 0.943 au. For 2I/Borisov, we quote values based on our CX emission model, but adopted 
the measured Q(CO) (Bodewits et al. 2020). For both objects we also estimate Lx from an empirically derived relationship 
with Lopt- 


Object Method R (au)  Qga (s!)  Lx(ergss |) Fx (ergss ! cm ?) 
2I/Borisov Measured Q(CO) and modeled CX emission. 2 1.07 x 10?" 7.5 x 10"? 6.2 x 1071" 
2I/Borisov Empirical relationship between Lopt and Lx. 2 - 6.1 x 10? 6.0 x 10718 
1lI/'Oumuamua Modeled Q(CO) and CX emission. 0.943 — 2.0 x 10?" 1.2 x 10? 5.4 x 10 * 
lI//Oumuamua Modeled Q(CO) and CX emission at rperi. 0.25 2.8 x 10?5 3.3 x 107° 7.7 x 107" 
11/‘Oumuamua Modeled Q(N2) and CX emission. 0.943 1.3 x 1078 4.7 x 10 22x10 7? 
1lI/'Oumuamua Modeled Q(H2) and CX emission. 0.943 1.5 x 1075 6.9 x 104 3.2x10 7? 
1lI//Oumuamua Empirical relationship between Lop: and Lx. 0.943 - 2.6 x 10° 3.6 x 10717 


was anomalous with respect to other orbital parameters, 
such as distance from the solar apex at encounter, and 
eccentricity). By contrast, 1I1/‘Oumuamua’s perihelion 
distance of 0.25 au was unexpectedly small. 

Next, we sample Earth’s position as a unit vector with 
azimuthal and polar angles jointly distributed according 
to p(0, 9) = sin ġ/4r, and scaled by 1 au. In reality, the 
orientation of the ISO's trajectory with respect to the 
ecliptic is determined by the inbound velocity vector and 
impact parameter. Stellar kinematics are a good proxy 
for the galactic distribution of ISOs; however, the typical 
ages of ISOs are poorly constrained at a population level. 


We assume the above distribution of (0, ¢) for our order- 
of-magnitude calculations. The distances R* = ||R(¢*)|| 
and Až = ||R(t*)—R’(t*)|| are solved for numerically by 
fixing Earth's position and assuming ~2"" = 70°. These 
distances are shown in Figure 3, with 1I/‘Oumuamua’s 
trajectory as an example. In cases where the entirety of 
the ISO's orbit lies within the solar avoidance zone, the 
Monte Carlo X-ray flux is immediately set to zero. 

We neglect the hyperbolic velocity vo; in our anal- 
ysis since the vast majority of ISOs accessible to the 
Rubin Observatory will have vo < 40 km s71. Larger 
velocities, however, would augment the solar wind ion 


collision rate by a factor of 10% or greater. It is prudent 
to conduct X-ray observations of such ISOs during the 
inbound portion of their trajectories. 


4.2. Bulk Properties 


The composition of a neutral medium does not 
strongly affect the X-ray power density from CX in the 
optically thick regime. However, the coma's density pro- 
file depends on the neutral production rate, which in 
turn depends on the volatility of surface material in the 
nucleus of an ISO. While Equation B11 estimates neu- 
tral production for bulk compositions similar to C/1996 
B2 (Hyakutake)'s, it is straightforward to account for al- 
ternative materials. The flux of molecules ejected from 
a unit surface area patch of the object (in units of s^! 
cm ?) is (e.g., Seligman et al. 2022), 


_ 1 - p) = eoT Sy, " 

AH/NA t ykp Tsu 
In the above, p is the albedo, I(t) is the time-dependent, 
solar irradiance, and e is the surface emissivity. Besides 
I(t) and physical constants, each parameter in the above 
equation is material-dependent. In practice, grey-body 
emission from the minor body is negligible and the in- 
cident radiation term dominates. Using the definition 
in Equation B12, exchanging I(t) for the incident solar 
flux, and integrating Equation 9 over the illuminated 
surface area (Xz) yields, 


Qs -ü-nf[ (pq) er w 


Albedo is unconstrained but is likely of order p ~ 0.1 
for most minor bodies in the Solar System. The av- 
erage projected surface area over all viewing angles is 
1/4 times the total surface area for any convex object 
(Meltzer 1949), so X; ~ X/4, where X is the total sur- 
face area of the object. For a sphere, Y; = 1a”, where 
a is the radius as defined earlier in this paper. There- 
fore, Equation 10 has the same functional form as Equa- 
tion B11. 

Hoover et al. (2022) noted that the size distribution of 
ISOs could enhance the Rubin Observatory’s yield and 
estimated detection rates for a range of absolute magni- 
tudes (see their Table 2). Many minor body populations 
adhere to a power-law size distribution 


N 


n(> a) = ka, (11) 


for some powerlaw index 8 and normalization constant 
k. Do et al. (2018) note that 6 > 3 for ISOs, or else 
the mass would diverge at large radii. In a collision- 
ally evolved system 8 = 2.5 (Dohnanyi 1969). However, 
Solar System comets have a more complex collisional 
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history, with a powerlaw index that depends on the size 
regime. For example, 8 = 1.45 in the 1 — 10 km range, 
and 6 = 1.91 in the 2— 5 km subset (Meech et al. 2004). 
For main-belt asteroids, 6 = 1.0 — 1.4 for sub-km di- 
ameters (depending on location within the belt), and 
B = 1.8 for larger objects (Yoshida et al. 2003). The 
distribution for ISOs may be shallow in the case that 
they are ejected from their respective systems before at- 
taining collisional equilibrium. Also, the bulk of ISOs 
may be well-described by a single value of 6, but the 
frequency will necessarily drop off for the largest and 
smallest objects due to physical limitations. We incor- 
porate the uncertainty surrounding 6 into our analy- 
sis by drawing it randomly from a uniform distribution: 
B ~ U(1.0, 2.0). As more ISOs are detected, it will be 
important to reconsider their size distribution. 


4.3. X-ray Flur 


We further restrict the number of free parameters by 
assuming the Rubin Observatory will detect 1.5 ISOs 
per year, the average of the optimistic and conservative 
rates presented by Hoover et al. (2022). The remaining 
free parameters are a, R, and composition. We ran- 
domly draw a nuclear radius a by first drawing a value 
for B, and subsequently sampling Equation 11. We trun- 
cate the distribution at agi, = 25 m and amax = 10 
km, which prohibits ISOs larger than most Solar Sys- 
tem comets and also guarantees they are large enough to 
realistically be detected by the Rubin Observatory (an 
object with 1/4 the radius of 1I/'Oumuamua and with 
the same albedo will have V — 24 at the same q — 0.25 
au). The new cumulative probability distribution is: 
(a-f nu) 


max 


MOOS Baan) P 
The minimum allowed heliocentric distance R = R* is 
determined per Equation 7, and the neutral production 
is calculated using Equation B11 for four representative 
pure compositions: H5, Nə, CO, and H20. We make an 
important modification, however, by drawing the helio- 
centric distance scaling exponent from a random uni- 
form distribution x ~ U(—4.33, —1.82) as opposed to a 
fixed x — —2, per Combi et al. (2019). The CX flux 
is then calculated using Equation B6 and an optically 
thick cross-section given by Equation B9 and correction 
factor &. 

Figure 5 summarizes our results, and demonstrates 
the viability of using X-rays as a diagnostic probe of 
ISOs across different compositional assumptions. As 
discussed earlier, composition has a profound effect on 
the extent of the outgassed coma and the expected 
X-ray flux. In the most conservative scenario that 
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ISOs are predominantly composed of water ice, then 
about 3% found by the Rubin Observatory will exhibit 
Fx = 2.7 x 1071? erg s^! cm~?, which can be robustly 
detected by a 10 ks exposure (ie., Tss = 10 ks cor- 
responds to a 5o detection; about 6% are accessible 
with a 75, = 100 ks exposure). Another 90% of H20- 
dominated ISOs will exhibit undetectable levels of CX 
emission (ie., Tss 21000 ks). If CO is the dominant 
constituent of ISOs, then about 1096 will emit X-rays 
detectable with a 10 ks exposure. The detectability 
fraction is even higher for the most volatile composi- 
tions, reaching 31% for Nə ice and 34% for Hə ice. With 
a longer 100 ksec exposure, these fractions reach 1996, 
47%, and 50% for CO, Ng, and Hs, respectively. 

These statistics suggest that the Rubin Observatory 
will discover new ISOs amenable to X-ray spectroscopic 
characterization. If the survey finds 15 2I/Borisov ana- 
logues (i.e., CO-dominated interstellar comets) over a 
ten-year campaign, then approximately 1 — 2 of them 
will be accessible with a 10 ks EPIC pn exposure. The 
Rubin Observatory will probably detect one such ISO 
in its first five years of operation. The approach is 
similarly effective for 1I/‘Oumuamua analogues. X-ray 
emission should be detectable for about one third of Nə- 
dominated and H5-dominated ISOs. On the other hand, 
consistent non-detections of CX emission will point to an 
alternative, refractory composition. For objects that ex- 
hibit nongravitational acceleration, lack of X-rays may 
also indicate an alternative mechanism that provides the 
required force (e.g., radiation pressure). 


5. CONCLUSIONS 


In this paper, we propose X-ray spectroscopy as 
a means of detecting outgassed comae and revealing 
highly volatile compositions of interstellar objects, a new 
class of astrophysical phenomena. The neutral gas co- 
mae of ISOs should undergo charge exchange with solar 
wind ions in the same manner as has been widely ob- 
served with Solar System comets. We suggest a tiered 
list of scientific goals that are realistically achievable by 
existing and concept X-ray missions, sorted by feasibil- 
ity: 


1. Photometric detection of CX (ie., at > 5o sig- 
nificance) in order to detect an outgassed coma. 
Modest exposure times of order 10 ks may reveal a 
statistically significant X-ray flux originating from 
an interstellar object, particularly when observed 
at heliocentric distances R < 1 au. The viable 
parameter space is depicted in Figure 6 for pn on- 
board XMM-Newton. Confidently detecting X-ray 
emission over a bandpass of 0.5 — 0.7 keV would 
immediately reveal the existence of a coma under- 
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Figure 5. Cumulative distribution function of X-ray flux 
(Fx) for ISOs detected by the Rubin Observatory. Equiva- 
lently, the fraction of ISOs detected by the Rubin Observa- 
tory that will exhibit a given X-ray flux or less. Four distri- 
butions are depicted, each corresponding to an ISO popula- 
tion comprised of a different pure volatile: H2O (blue), CO 
(red), N2 (green), or H2 (orange). On the top x-axis, Fx is 
converted into an exposure duration (Tso) that would yield 
a 5o detection with XMM-Newton EPIC pn. The vertical 
blue line delineates Fx c 2.7 x 10^ !? erg s^! cm ?; a sig- 
nal which would be detected at 5o confidence with at 10 ks 
exposure. For example, about 1096 of CO-dominated ISOs 
detected by the Rubin Observatory will exhibit this X-ray 
flux or greater. 


going CX with solar wind ions. X-ray observa- 
tions of 1I/‘Oumuamua near its close approach to 
Earth would have either confirmed the presence 
of a coma, or placed an upper bound on the pro- 
duction rate to a point that another mechanism 
would be necessary to explain the anomalous ac- 
celeration. 


2. Joint constraints on a and Qgas in order to infer an 
ISO's volatile content. Once an ISO's ephemeris 
is determined, the X-ray flux is degenerate with a 
and H. By assuming a range of plausible albe- 
dos, one can constrain a from an optical light 
curve (e.g., Mashchenko 2019). The radius is also 
constrained by the energetics of nongravitational 
acceleration. For example, Seligman & Laughlin 
(2020) suggested various ices (Ho, Ne, No, Ar, 


O2, Kr, Xe, CO2, H20) that could be compati- 
ble with 1I/‘Oumuamua’s outbound acceleration. 
For a subset of these species (No, H2, CO, Ne, 
plus CH4), Jackson & Desch (2021) showed the re- 
lationship between nongravitational acceleration, 
albedo, and mean spherical radius. Across these 
different species, the allowable range for a lies be- 
tween 20— 25 m for albedo p > 0.6. However, their 
respective sublimation energies (H) vary by a fac- 
tor of 10, which, according to our model, implies 
A varies by a factor of 100. This interaction area 
is probed by the X-ray flux from CX. If the nu- 
clear radius of an ISO can be constrained within 
a factor of order unity, then an anomalously high 
Fx would immediately reveal an exotic, highly- 
volatile composition due to the inverse proportion- 
ality between Qgas and H. We reiterate that a 
non-detection in X-ray would favor low-volatility 
ices (e.g., H20, CO2) or an alternative accelera- 
tion mechanism. 


3. Measuring line strength ratios towards robustly 
determining the coma composition. Mullen et al. 
(2017) demonstrated that relative CX emission 
line strengths (including strong transitions from 
Cv, Cvr, and Ovi) depend on the neutral tar- 
get medium. Specifically, No, H20, O, OH, CO, 
and CO» targets were considered in their study. A 
high-fidelity X-ray spectrum may be achievable by 
existing facilities for the brightest ISOs. However, 
this method is most amenable to next-generation 
X-ray telescopes which feature both large effec- 
tive areas and high-spectral resolution. A detailed 
study of this approach, and its feasibility with X- 
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ray concept missions, would be scientifically use- 
ful. 


Baseline estimates of the Rubin Observatory's yield 
suggest around 15 new ISOs over a ten-year campaign. If 
these objects are predominantly composed of CO, about 
one or two will exhibit detectable levels of X-ray flux. 
In the case that most ISOs are comprised of Hə ice, CX 
will be observable for about five of them. Efforts to un- 
derstand this mysterious population benefit from broad- 
band spectroscopy — infrared, optical, and UV observa- 
tions provide complementary insights into physical and 
chemical properties of the nucleus and coma. As X-ray 
observations of CX with comets have been an active area 
of study, similar efforts directed toward ISOs will at the 
very least reaffirm findings made at other wavelengths, 
and potentially discern the true compositions of these 
minor bodies. 
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APPENDIX 


A. REVIEW OF SOLAR WIND CHARGE EXCHANGE 


Since the discovery of strong X-ray emission from a cool cometary source, C/Hyakutake 1996B2, in 1996 by Lisse et 
al., solar wind charge exchange has gained considerable traction for its ability to explain X-ray emission from objects 
within the Solar System. CX models have successfully explained observed X-rays from the volatile comae of Solar 
System comets (Cravens 1997; Wegmann et al. 1998), as well as emission originating from Earth (Wargelin et al. 2004; 
Fujimoto et al. 2007), Mars (Kallio et al. 1997; Dennerl 2002), Venus (Dennerl et al. 2002), and Jupiter (Metzger et al. 


1983; Cravens et al. 2003). 


In order to produce X-ray emission, heavy ions from the solar wind interact with cool atomic and molecular gas 
surrounding the object. This process also produces a source of X-ray emission from star-forming galaxies such as 
M82 (e.g., Liu et al. 2012; Zhang et al. 2014). In such a galaxy, CX occurs at the interface between hot and cool 
interstellar gases, likely responsible for the enhanced forbidden transitions of Ka triplet emissions from He-like ions 
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and high-order transitions in Lyman series of H-like ions (Zhang et al. 2014). There has also been tentative detection 
of charge exchange in galaxy clusters (Gu et al. 2018). 
The general CX reaction is 
AIt +N > CAPSISHDHE 4 NF, (A1) 


where the ion A9* gains an electron from the neutral species N (atomic or molecular, often H, H5 and He). The integer 
q denotes the initial ionization state of the species A. The electron is in an excited state, denoted by the superscript 
*. As it cascades to the ground state, the electron emits X-ray and UV photons. Implicit in this equation are initial 
and final state quantum numbers (n, l, and m) of the transferred electron. The most probable final-state principal 
quantum number n is given by Equation 2.6 of Janev & Winter (1985), which is an approximation that is generally 
applied while modeling this process (Smith et al. 2012). 

The CX interaction cross-section (054) for a given species (denoted by subscript s) and charge (denoted by subscript 
q) was modeled by Wegmann et al. (1998), and is given by, 


q-1 
Osq = 
q q?/2n? — |15,,/27.2eV | 


2 
) x 0.88 x 1016 (cm?). (A2) 


The variables q and Jp, in Equation A2 correspond to the dimensionless integer charge of the ion and ionization 
potential (measured in eV) of the target species, respectively. 

For comets, the solar wind typically interacts with volatile coma particles. Some of the most common target species 
outgassed by Solar System comets are water and its constituents (O, OH, H) (Cochran et al. 2015; Biver & Bockelée- 
Morvan 2016; Bockelée-Morvan & Biver 2017). All of these species have similar ionization energies ~ 13 eV (Wegmann 
et al. 1998). CX cross-sections are ~ 107! to 10714 cm? for heavy solar wind ions such as Ovi, Ovi, Cvi, Cvr, 
Nvu, Nerx, Six, and Fexir. That is, charge exchange occurs when these ions come within ~ 1 nm of the neutral 
targets. The energy difference between the captured electron's initial (excited) state and ground state (Fexcit) can 
reach several hundred eV. Therefore, these transitions provide significant contributions to EUV and X-ray radiation. 
CX emission from cometary comae has rich observational and theoretical backing, as reviewed by Bodewits et al. 
(2012) and Dennerl et al. (2012). 

Our main concern is charge exchange between the solar wind and the outgassed species within an ISO's coma. 
Assuming that a given solar wind ion undergoes charge exchange once during its passage through the neutral medium, 
the X-ray power density P4 for transition (j) at a given position r is: 


Ps (r) = Psqlr) Osq bsaj Na (r) AEsaj , (A3) 


(Cravens 1997, 2000; Cravens et al. 2009) which has units of eV cm ? s^!. Here ®,, denotes the solar wind flux 
of a given ion, sq is the CX cross-section (Equation A2), bsqj is the spectral cascading probability for transition (j) 
which releases a photon of energy AFggj, and n, is the number density of neutral species. The solar wind flux may 
be written as $44 = fiqnswusw, where ngw and usw correspond to solar wind proton number density and velocity, 
respectively. The solar wind fraction of a given species/charge is fzq. These quantities are functions of position in 
space. For example, the unshocked solar wind has a typical speed of ~ 400 km s^! and density of ~ 0.4 cm^? at 5 au 
(Cravens et al. 2003). The density nsw falls off as the square of distance from the host star (Cravens 2000). While 
the solar wind flux varies over time, the fraction of heavy ions to protons in the solar wind remains fairly consistent at 
f 7 107°. Intensity is obtained by integrating Equation A3 over a given path length (i.e., performing a line integral). 


B. CHARGE EXCHANGE WITH OUTGASSED COMAE 


Our prescription for CX emission is calibrated to observations of C/1996 B2 (Hyakutake) and early models that 
described its X-ray luminosity. Wegmann et al. (1998) calculated the total X-ray emissivity by summing over all solar 
wind ions in Equation A3. They adopted parameter values of a4, = 3 x 1071? em?, fo = 0.0005, and 1100 eV emitted 
per oxygen ion in the solar wind based on summing the CX excited state energies of all ions and weighting them 
according to their relative abundance fraction. Further, they assumed an ‘effectivity’ of 0.4 which replaces the b; 
terms. The power density is approximately 


Px(r) = 4 x 10 7? nswn,(r')ergem?s ! , (B4) 


where the density of neutrals at distance r’ from the comet’s nucleus follows 
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TS ( Qes ) S (B5) 
Arv|r’| 

for a total neutral production rate Qeas, outflow velocity v, and photodestruction length scale A. Wegmann et al. 

(1998) adopted Qgas = 1.5 x 107? s^! and nsw = 7 cm ?. By neglecting photodestruction, which Cravens (1997) 

deemed negligible within 5 x 10° km, and adopting a nominal v = 1 km s71, integration over a sphere of beam radius 

135, 000 km yields Lx ~ 6 x 1015 erg s~!, comparable to the measured luminosity and other estimates in the literature. 

As follows, we determine expectations for CX emission from the comae of other minor bodies. 


B.1. Luminosity from the Optically Thick Coma 


CX produces a surface luminosity, 4x1 (Cravens et al. 2003), given by the equation, 


4nI = nsw usw f NAE. (B6) 


Each ion contributes, on average, N photons of typical energy AE, and the solar wind minor ion fraction is given by 
f. The total luminosity is obtained by integrating Equation B6 over the effective surface area. Our baseline model 
assumes that this area is the region where the coma is optically thick to CX. That is, it satisfies the criterion: 


| ( vaga (1^) )da> 1 (B7) 
S 


for a linear projected path S of the solar wind through the coma, where r' is the distance from the cometary nucleus 
as in Equation B5. As long as the coma neutral column is large enough that the solar wind is fully depleted of all 
highly stripped minor ions, the individual transition probabilities are not important. For simplicity, we also fix the 
cross-section at Osą = 3 x 1071 cm?. This value is consistent with expectations from the ionization potentials of Hə 
(Ip,s © 15.4 eV) and No (Ip, © 15.6 eV). Considering these compositions and those more typical of comets, we find 
2x 1071? cm? € o4 < 13 x 10? cm? for ion charges +5 € q X +7. We assume that the coma is approximately 
spherical and follows a neutral density profile given by Equation B5. We define the impact parameter, b, as the 
minimum distance between the path S and the nucleus. Also, let R be the distance between the Sun and the nucleus, 
and assume that R >> b. We parameterize the path S as a function of r’. We define the dimensionless parameter 
x = |r’|/b and a coefficient K = osqQgas/4nw. By substituting Equation B5 into Equation B7 with these new 
parameters, the optically thick transition occurs where 


2K [? 1 hf 
e "dy =1, B8 
b J zyr? — 1 G 
or, equivalently, where 
T/2 
=| eo sec8/A ag E (B9) 
b Jo 


Equation B9 is equivalent to a function b = b(a,4, W, A, Qgas) Which may be evaluated numerically. It returns a distance 
within which the coma is optically thick to CX. For physically plausible photodestruction scales of A = 105 — 10!! cm 
(Combi et al. 2004), a reasonable approximation is that b ox Qgas in the regime Qgas < 1078 s71. Figure 6 shows the 
validity of this approximation for different assumed A. As A — oo, the function approaches a linear relationship. 

Taking C/1996 B2 (Hyakutake) as a nominal example (hence the above assumed A and Qgas values), then we find 
b = 10,000 km which is in reasonable agreement with the 30,000 km penetration depth estimated by Wegmann et al. 
(1998). In order to estimate total X-ray luminosity, we assume nsw = 7 cm? at 1 au, usw = 400 km s^!, f = 0.001, 
N = 1, and AE = 550 eV. The above equations predict Lx = 8.2 x 10!^ erg s^! from the optically thick region 
alone. This estimate is about 5x less than the total luminosity measured by ROSAT. Using a refined hydrodynamical 
model, Wegmann et al. (2004) also estimated Lx for C/1996 B2 (Hyakutake) during its close approach to Earth. Their 
measurements were centered around 101° erg s~!, which is higher than the initial estimate from Lisse et al. (1996). 
The discrepancy between these literature estimates and our own is probably due to significant emission that took place 
in optically thin regions of the coma, currently unaccounted for in our model for the optically thick zone. Therefore, 
it is useful to calibrate our model and correct for this component per the following. 
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Figure 6. Relationship between outgassing rate Qgas and radius b within which a coma is optically thick to CX. Since b lacks 
a closed form (Equation B9), it is solved for numerically. The function is plotted for several assumptions of the ionization scale 
A. In all cases v = 1 km s^!. Our main analysis assumes À = 5 x 101° cm (red line). As A increases, the function approaches a 
linear relationship between Qgas and b (black dashed line). 


B.2. Extensions to the Analytic Model 


Holding other parameters constant, b scales linearly with Qgas in the regime Qgas < 10°° s^! for A ~ 5 x 101° cm. 
The CX interaction area in this regime can be approximated by 7b? zz TQ%,,, with P = 1.5 x 1077? s? cm?. A similar 
conclusion was reached by Wegmann et al. (2004), who developed a hydrodynamical model of X-ray emission morphol- 
ogy and applied it to observations of comets. They derived a relationship Lx = CH oe where H = uswnsw f NAE 
is the heavy ion solar wind flux. By fitting their estimates of H os to observed X-ray luminosities of comets, they 
determined a constant of proportionality C = 10738 s? cm?. Let A define an effective, optically thick surface area 
that would yield a given Lx. Then the relationship found by Wegmann et al. (2004) is equivalent to A = C a We 


adopt a correction factor k = y C/T ~ 8.1 that defines a new effective radius, 


b 


Kb (B10) 


In other words, & is a factor that increases the radius of the CX interaction area from b to b, thereby approximately 
accounting for emission from regions beyond the optically thick zone. Prior observations of CX within cometary comae 
(e.g. Lisse et al. 1999; Lisse et al. 2004) have produced detections of significant X-ray emission extending as little as 
~ 5 x 10* km (for 2P /Encke 2003), and as far as an apparent limit at ~ 10° km (for C/1996 B2 (Hyakutake), C/1991 
K1 (Levy), and C/1995 O1 (Hale-Bopp)). 

We adopt Qeas c a?R ?, where a is body's radius, and R is the instantaneous heliocentric distance. The R^? 
dependency is the same as that of the incident solar radiation. While physically motivated, this scaling simplifies 
true cometary activity. Following a long-baseline survey of water outgassed by 61 comets, Combi et al. (2019) found 
generally steeper scaling: Q(H20) x R”, with —4.33 € x < —1.82 depending on taxonomic class. This finding is 
incorporated into our expectations for future ISO observations in 84. 

Our model must also account for the fact that more volatile species will sublimate and outgas at higher rates 
(potential compositions are discussed further in §4). We use the following scaling relation, 


Qs = 1.5 x 10% x a) Gu) (= " icm a) 5^ (B11) 
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which is calibrated to the nuclear radius of C/1996 B2 (Hyakutake) (Lisse et al. 1999), and the production rate 
adopted by Wegmann et al. (1998). Production rates for OH were measured by Schleicher (1996), and varied from 
log,g Q(OH) = 28.89 (at R = 0.94 au) to log;, Q(OH) = 29.17 (at R = 1.08 au). In the above, H is the total energy 
input for each coma particle, 


H = AH/Na + ykTsup - (B12) 


The above expression involves enthalpy of sublimation (AH), temperature of sublimation (Tsu), and adiabatic index 
of the escaping vapor (^), all of which are material dependent (Table 1). The quantity H is the sum of two components: 
the energy required to sublimate an ice molecule (the first term) and the energy required to heat the molecule to the 
outflow velocity (the second term). In the case of water, Hy,0 = 9.3 x 1077? kJ (values for other species are listed in 
Table 1). Chemical properties also determine the outgassing velocity v: 


vc e, =  ykTsay/umg , (B13) 


for mean molecular weight u and sound speed cs. For typical values v = 0.3—1 km s71, Q < 108° s~!, photodestruction 
timescales tphoto < 107 s, and A= wtphoto, the optically thick radius does not exceed 10° km under our model. It is 
limited by the exponential decay term in Equation B5. We reiterate that only monotomic and homonuclear diatomic 
species (having zero diplole moment) will be more amenable to CX X-ray characterization than to infrared fluorescence 
spectroscopy. 

The normalization in Equation B11 is to C/1996 B2 (Hyakutake), and the expression should be considered an 
order-of-magnitude estimate of Qeas for other objects. Ideally, one may use independent estimates of Q4; obtained 
either from infrared spectroscopy or a more detailed theoretical model that accounts for the target's composition. For 
most Solar System comets, Qeas ^ Q(H20) + Q(CO2) + Q(CO), which is the sum of production rates of the most 
common volatile components in comets (Bockelée-Morvan & Biver 2017). Nevertheless, Equation B11 enables baseline 
predictions for new minor bodies, and is especially useful for population-level predictions when given distributions for 
a and perihelion rye. Additionally, we assume that the solar wind proton number density is given by, 


ngw(R) = 7 cm? I (B14) 


which is maximized at perihelion when R = rpe,j. Note, flares can enhance the solar wind ion flux, thus improving the 
feasibility of CX X-ray observations at larger R. For example, a flare nearly doubled the soft X-ray count rate from 
C/1999 S4 (LINEAR) (Lisse et al. 2001). 

Combining Equations B6 & B10, we arrive at a general model for X-ray luminosity 


Lx = rb ngwusw f NAE, (B15) 


which is a function of {0sq, A, a, rperi, 4, AH, y, Tsun}, assuming observations take place at perihelion. Also, a nominal 
A = 5 x 101° cm is adopted for our analysis. We use this model to predict the luminosity of CX emission with ISOs 
in 84, and a schematic of the entire CX process with ISOs is shown in Figure 1. It is worth explicitly highlighting 
our model's dependence on perihelion distance. We demonstrated an approximately linear relationship between b and 


Qeas, where Qaa; cx Eme through its dependence on the solar radiation flux Fo « ms 'The CX surface area follows 


A x b?, and the incident heavy ion flux follows nsw ox prom Therefore, Lx c (enim the perihelion of an ISO strongly 


dictates whether its CX X-ray emission is detectable (similar scaling relationships were explored by Lisse et al. 1999). 
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